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MERIDIONAL TILT OF THE STELLAR VELOCITY ELLIPSOID DURING BAR BUCKLING INSTABILITY 
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ABSTRACT 

The structure and evolution of the stellar velocity ellipsoid plays an important role in shaping 
galaxies undergoing bar driven secular evolution and the eventual formation of a boxy /peanut bulge 
such as present in the Milky Way. Using collisionless N-body simulations, we show that during the 
formation of such a boxy/peanut bulge, the meridional shear stress of stars, which can be measured 
by the meridional tilt of the velocity ellipsoid, reaches a characteristic peak in its time evolution. It is 
shown that the onset of a bar buckling instability is closely connected to the maximum meridional tilt 
of the stellar velocity ellipsoid. Our findings bring new insight to this complex gravitational instability 
of the bar which complements the buckling instability studies based on orbital models. We briefly 
discuss the observed diagnostics of the stellar velocity ellipsoid during such a phenomenon. 
Subject headings: galaxies: bulges - galaxies:kinematics and dynamics - galaxies: structure - 
galaxiesrevolution - Galaxy: disk, galaxies:halos, stellar dynamics 



1. INTRODUCTION 

Understanding the structure and dynamics of a galaxy 
crucially depends on the knowledge of the three dimen- 
sional stellar distribution function (DF), which is not a 
direct observable. The first few moments of the DF, e.g., 
density, mean velocity and the velocity dispersion ten- 
sor together can provide important clues regarding the 
dynamical state of the galaxy and the gravitational in - 
stabilities it might have undergone ()van der KruitiH999l) . 
Of particular interest is the buckling instability of a 
stellar bar in a disk galaxy and the subsequent for- 
mation of a boxy /peanu t bulge dCombes et al.1 fT"990b 
Pfenniger fc Friedlil ll991t iRaha et al.1 119911 : iPfennigeri 
19931 IAthanassoulall2005| ). A bar buckles under its own 
self-gravity when it becomes sufficiently strong, thereby 
bringing substantial changes in the velocity distribution 
of stars and the galactic potential. One possible way to 
quantify such a change is to study the structure and evo- 
lution of the stellar velocity ellipsoid during the buckling 
instability and to provide potential diagnostic observ- 
ables. In particular, how is the tilt of the velocity ellip- 
soid related to th e boxy /peanut bu lge such as presented 
in the Milky Way ([Dwek et al.lll995l ). This requires, how- 
ever, an unambiguous identification of the buckling event 
that a galaxy might be experiencing. However, the on- 
set of buckling instability is not clearly understood be- 
cause it is difficult to follow the orbits of stars subject 
to a rapidly changing gravitational potential during the 
buckling. During this transient phase the dynamics is 
strongly collective and an orbit decomposition can only 
be a partial description of the process. Nevertheless, nu- 
merous simulation studies marked this event by a de- 
crease in the bar strength or in th e ratio of vertical-to- 
radial velocity dispersion (a z /( r r ) l)Combes et al.l H990t 
iMartinez-Valpuesta et all I2006D , providing a gross un- 
derstanding of this event. Simulations show that often 
such demarcation is blurred and a more precise indica- 



tor of this event would be useful and complementary to 
the already existing ideas. It is worth re-investigating 
the buckling instability and the relation it might have 
with the orientation of the stellar velocity ellipsoid, in 
particular with the tilt angle. 

The shape and orientation of the stellar velocity ellip- 
soid are tightly connected to the symmet r y of the un- 
derly i ng galaxy potential (iLindbladl 119301 iLvnden-Belll 
[19621 : lAihendt fc Cuddefordlll99in . In a stationary, ax- 
isymmctric disk galaxy the stellar velocity ellipsoid in 
the galactic midplane is perfectly aligned with the galac- 
tocentric coordinate axes, in other words, all the off- 
diago nal elements of the veloci ty dispersion tensor are 
zero (|Binnev fc Tremainel fl987t) . Thus, measuring the 
off-diagonal components of the dispersion tensor in ob- 
servation may provide one with an inference about the 
presence of non-axisymmetric features in a galaxy. Away 
from the galactic midplane, the tilt of the velocity el- 
liposid might depend on the mass distribution of the 
galactic disc as well as the flattening of the dark matter 
halo. In the context of the Mil ky Way, the analysi s of the 
RAVE survey data release 2 (jZwitter et al 2008) shows 
that t he velocity ellipsoid is tilted towards the Galactic 
Plane ([Siebert et al.ll2008l ) and has been nicely demon- 
strated in a recent paper by iPasetto et al.l (|2012f) . How- 
ever, the measured tilt angles can not put a strong con- 
straint on the disc parameters and halo flattening due to 
large proper motion errors and small sample size in the 
RAVE DR2 (jSiebert et al.l[200l . 

On the other hand, non-axisymmetric structures such 
as bars, spiral arms in disk galaxies might play an im- 
portant role in accounting for the observed orientation 
of the stellar velocity e llipsoid. Numerical study by 
IVorobvov fc Theisl (|2008f ) shows that the vertex devi- 
ation of the velocity ellipsoid is globally correlated to 
the amplitude of the spir al arms. Using Hipparcos data 
and dynamical modeling, iDehnenl (|2000f ) has shown how 
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the Galactic Bar (iBlitz fc Spergell 119911: iBinnev et all 
119911 : iDwek et"aH 11995ft could have influenced the ve- 
locity distribution in the solar neighborhood. The ob- 
served low-velocity streams in the solar neighborhood 
are also thought to h ave arisen due to the Galactic Bar 
(jMinchev et al.ll2010( ). It would be useful to understand 
how the presence of a bar or spiral arms which are highly 
time dependent would change the orientation of the ve- 
locity ellipsoid. 

We use N-body simulations to follow the evolution of 
the stellar velocity ellipsoid in a galaxy which undergoes 
bar instability and forms a boxy bulge during the secular 
evolution in a self-consistent way. The buckling of the bar 
causes the morphological evolution of the disk, convert- 
ing its central parts into a boxy /peanut bulge. In order 
to gain further i nsight into the physics of bar bucklin g 
(|MerriflelJl99fllMartinez-Valpuesta fc Shlosmanll2004f ). 
we investigate the role of anisotropic stellar pressure and 
show that there is a characteristic signature in the way 
the stellar velocity ellipsoid evolves. The primary goal of 
this paper is to understand the buckling event of a bar 
which forms the boxy /peanut bulge and its relation with 
the tilt of the stellar velocity ellipsoid. 

The paper is organized as follows. In the next section, 
we outline the general concept of the stellar velocity el- 
lipsoid and the relevant quantities that we measure from 
our simulation. Section [3] briefly describes the galaxy 
models used for the present study and simulation. The 
disk evolution and boxy bulge formation is described in 
section |U The shear stress and its relation to bar buck- 
ling is shown in section [5] We discuss the tilt angle of 
the velocity ellipsoid in section 15.11 Finally, section [5] 
presents the discussion and conclusions from this work. 

2. STELLAR VELOCITY ELLIPSOID 

The components of the velocity dispersion tensor at a 
radial location r in the stellar disk are computed from 
the velocity comp onents of a group of star s using the 
following formula (jBinnev fc Tremaindll987l ): 



and 



T ij = ( v i v j)-( v i)( v j)> (!) 
Uj denote the velocities of a group of 



where Vi 

stars. i,j = r,^,z in a cylindrical coordinate system. 
Angular bracket denotes the averaging over a group of 
stars. Given the velocity dispersion tensor, the stress 
tensor of the stellar fluid can be written as 



-p(r)a 2 



(2) 



where p(r) is the local volume density of stars at a posi- 
tion r. It is convenient to think of the entire stress tensor 
as a sum of two different kinds of forces acting on a small 
differential imaginary surface (dS) between two adjacent 
volumes of stars, i.e., 

t=t„+t s , (3) 

where t„^ = — /5(r)cr^ is called the normal stress acting 
along the normal to dS and cr^ are the diagonal com- 
ponents of the above matrix. r s ,y = —p(r)afj,i ^ j is 
called the shear stress, acting along a direction perpen- 
dicular to the normal to dS, i.e., in the plane dS. In gen- 
eral, the shape of the velocity ellipsoid is determined by 
the normal stress, and the shear stress is responsible for 
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Fig. 1. — Initial circular velocity curves for the model RCG051A 
and RHG057. In both panels, red dashed line denotes the disk, 
blue dotted line the bulge and green dash-dot line the dark halo. 
Solid black line denotes the total circular velocity curve. The inner 
regions of RCG051A are disk dominated, while that of RHG057 
are dark matter dominated. 



TABLE 1 

Initial disk, halo and bulge parameters. Q is the Toomre 

STABILITY PARAMETER AT 2.5R d ; M h AND M b ARE THE MASSES OF 
HALO AND BULGE. 



Models 


Q 


M d 


T or b 


M h /M d 


M b /M d 






(xl0 10 M Q ) 


(Gyr) 






RCG051A 


1.21 


4.75 


0.301 


6.52 


0.05 


RHG097 


1.84 


3.11 


0.306 


7.88 


0.43 


RHG057 


2.98 


0.86 


0.326 


22.53 


0.18 



the orientation or deformation of the ellipsoid w.r.t. the 
galactocentric axes (e r ,e v ,e z ). The orientation of the 
velocity ellipsoid can be computed using the off-diagonal 
components of the velocity dispersion tensor. The merid- 
ional tilt of the velocity ellipsoid can be computed using 
the following relation: 



©tilt = ^ arctan 



2aL 



(4) 



We evaluate the shape of the velocity ellipsoid and the tilt 
angle in the inner region of the disk where the dynamics 
of stars is dominated by a bar and study their evolution 
as the bar enters into the non-linear regime where an 
analytic understanding is difficult. 

3. INITIAL GALAXY MODELS 

In order to study the evolution of the stellar velocity 
ellipsoid subject to a non-axisymmetric bar potential, we 
perform a large number of s imulations of isolated galax - 
ics built using the method of lKuiiken k, Dubinskil (|1995f) . 
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Of these, we present here 3 fiducial models (named as 
RCG051A, RHG057 and RHG097) of disk galaxies with 
varying dark matter distribution and Toomre stability 
parameter (Q). The initial disk has an exponentially de- 
clining surface density with a scale length R c i and mass 
Md- The live dark matter halo and bulge are modelled 
with a lowered Evans and King DF respectively. For fur- 
th er details on mode l const ruction, the reader is referred 
to ISaha et all (|2010l 120121) . We scale the models such 
that Rd = 4 kpc and the disk masses are given in Ta- 
bled! Orbital time scales T or b (at 2.5Rd) and other initial 
parameters are given in Table [1] In Fig. [TJ we show the 
circular velocity curves for RCG051A and RHG057. Cir- 
cu lar velocity curve for RHG097 can be found in Fig. 2 
of lSaha et all (|2010h . 

The simulations we re performed using the Gadget code 
(jSpringel et al.ll2001[ ) which uses a variant of the leapfrog 
method for the time integration. The gravitational forces 
between the particles arc calculated using the Barnes- 
Hut tree algorithm with a tolerance parameter t oi = 
0.7. The integration time step used was ~ 0.82 Myr 
for RCG051A, 0.65 Myr for RHG097 and 1.5 Myr for 
RHG057. Two of these models were evolved for about 
6 — 7 Gyr, and RHG057 was evolved for about 12 Gyr to 
understand the long term evolution, bar growth and the 
asymptotic properties of the stellar velocity ellipsoid. 

Each of these models were constructed using a total 
of 2.2 million particles, out of which disk and halo have 
1.05 million each and 0.1 million particles are assigned 
to the bulge. The softening lengths for disk, bulge and 
halo particles were chosen so that the maximum force 
on ea ch particle is nearly the same (jMcMillan fe Dehnenl 
120071) . In the model RHG097, the softening lengths used 
for the disk, bulge and halo were 12, 25 and 33 pc respec- 
tively. For RCG051A, they were 12, 10 and 31 pc and 
for RHG057, 12, 17 and 57 pc respectively. The total 
energy is conserved well within 0.2% till the end of the 
simulation. 

4. DISK EVOLUTION THROUGH BAR GROWTH 

Bar driven secular evolution is an important internal 
process through which galaxies change their morphol- 
ogy and kinematics. The rapidity of such a process de- 
pends on various factors of which bar strength plays a 
significant role. A bar forms out of the disk instabil- 
ity and grows via nonlinear processes as the disk stars 
exchange energy and angular momentum with the sur- 
roundi ng dark matter h alo and a preexisting classical 
bulge (|Saha et al.l l2012f ) through gravitational interac- 
tion. The evolution of an initially axisymmctric stellar 
disk and growth of a bar is depicted in Fig. [2j Fig. [3] and 
Fig. |4] which present the surface density maps of all the 
stars including that of a preexisting classical bulge. 

In Fig.[5j we present the time evolution of the bar am- 
plitude measured by the m = 2 Fourier component of the 
surface density of disk stars alone for the three fiducial 
models mentioned above. The growth rates of bars are 
significantly different in these models which also differ in 
the relative fraction of dark matter within the disk re- 
gion. In model RHG057, the dark matter dominates the 
disk right from the center of t he galaxy, see Fig . [TJ Ac- 
cording to the classification of ISaha et al.l (|2010f ). model 
RHG057 forms a type-II bar and models RCG051A and 
RHG097 form a type-I bar. Typically, type-I bars are 
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Fig. 2. — Face-on surface density maps of all the stars belonging 
to the disk and a preexisting classical bulge in the galaxy model 
RCG051A. Top left panel shows surface density at T = 0, top right 
at 2.0, bottom left at 3.4, and bottom right at 5.5 Gyrs. 
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Fig. 3.— Same as in Fig.[2]but for RHG097. Top left panel shows 
surface density at T = 0, top right at 2.0, bottom left at 3.4, and 
bottom right at 5.5 Gyrs. 

strong and go through the well known vertical buck- 
ling instability (|Combes et al.|[l9 90; Pfe nniger fc Friedlil 
19911: iRaha et al.l I199U iMartinez-Valpuesta et al.l 120061: 
Debattista et alj ' 20061) leading to the formation of a 
boxy/peanut (hereafter b/p) bulge as depicted in Fig.lTTI 
and Fig. [12] Whereas type-II bars which are weak and 
grow on secular evolution time scale, normally do not 
go through any appreciable buckling instability. We 
evolved the model RHG057 for a Hubble time and the 
disk showed no signature of buckling instability, although 
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Fig. 4. — Same as in Fig.[2]but for RHG057. Top left panel shows 
surface density at T = 0, top right at 2.0, bottom left at 3.5, and 
bottom right at 7 Gyrs. 
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Fig. 5. — Time evolution of normalized bar amplitude for model 
RCG051A (green), RHG097 (red) and RHG057 (blue). 

it has grown a moderate size bar by that time. As a re- 
sult, the disk in this model did not form a b/p bulge 
by that time. A thorough understanding of the buckling 
instability would perhaps require a tool combining the 
orbital analysis and collective effect of the stars in the 
disk and their role at the 2 : 1 vertical inner Lindblad 
resonance (ILR). In order to pinpoint the location of the 
ILR, corotation resonances (CR), we first compute the 
disk frequencies (0, k, v) by a direct sum of the first and 
second derivatives of the A-body potential obtained from 
the reflection symmetrized particle distribution with re- 
spect to the z = plane and the R = rotation axis. 
We carry out this exercise for each snapshot and com- 
pare n and fi — v/2 with the pattern speed of the bar. 
Fig. |5] shows the locations of the vertical ILRs and CRs 
at two different epochs for the two models RHG097 and 
RHG057. In the case of RHG097, the locations of verti- 
cal ILRs before (t = 2.2 Gyr) and after (t = 4.8 Gyr) the 
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Fig. 6. — Vertical resonances in the stellar disks of two galaxy 
models mentioned in the figure. Green dashed lines represent the 
angular frequencies (Q) at t = 2.2 Gyr (for upper panel) and t = 
3.6 Gyr (lower panel). In the upper panel, Q — u/2 profiles (red 
solid and blue dash-dot lines) are drawn at t = 2.2 and 4.8 Gyr and 
the corresponding bar pattern speeds at those times are denoted 
by the dashed and solid black lines. In the lower panel, they are 
at t = 3.6 and 7.0 Gyr and the corresponding bar pattern speeds 
are denoted by horizontal black lines. The unit of frequencies for 
RHG097 and RHG057 arc 44.65 and 19.0 kms -1 kpc -1 . 

buckling instability arc still within the bar region indi- 
cating that the orbits lie close to the 2 : 1 vertical oscilla- 
tions before and a fter the peanut formation. A detailed 
orbital analysis bv lPfenniger fc Friedlil (|1991|) shows that 
the 2 : 1 vertical resonance is essentially responsible for 
the formation of b/p bulge, but the collective behaviour 
of stars in the vicinity of such resonances remains ob- 
scured. For example, it is not understood what is the 
role of shear stress or the anisotropic stellar pressure in 
such a process which eventually leads to the formation of 
b/p bulges. Below we elaborate on the possible relation 
between the buckling instability and the structure and 
evolution of the stellar velocity ellipsoid, in particular 
the shear stress due to the disk stars. 

5. SHEAR STRESS AND BUCKLING INSTABILITY 

As the bar becomes stronger it enters into the regime 
of the buckling instability. This instability is highly 
nonlinear in the sense that the bending oscillation am- 
plitude is not proportional to self-gravity. In the case 
of a bending instability, we would expect a propor- 
tionality between the bending oscillation and the im- 
posed load (here, self- gravity). In the linear regime, 
a stellar disk is stable against the low order (e.g., 
rn = 0, 1) bending perturbation and the self-gravity of 
the perturbation acts like a stabilizing a gent as show n 
by several authors, e.g., Toomre (1966 1 ). lArakil (|1985l) . 
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Fig. 7. — Time evolution of A\^ z denoting the vertical asymmetry 
in the bar region. Green line denotes model RCG051A, red line for 
RHG097 and blue for RHG057. 



IMerritt fe Sellwoodl (fl99l . lSellwoodl (fl996h . lSaha fe Jod 
(|2006l) . It is the pressure forces which destabilizes 
a stellar disk in response to a ben ding p e rturba tion. 
From the works of Toomre ( 1966), lArakil (|1985t l and 
iFridman fe Pol iachcnkol (|1984[ ). we learnt that a stel- 
lar slab of finite thickness would go bending unstable 
if <J z ja r < 0.3. However, the critical value of a z /a r , 
at which a self-consistent rotating stellar bar would go 
bending unstable is unclear. Actually, it is doubtful that 
a criterion based only on a local quantity such as <J z /a r 
would apply in a bar, as the 2/1 vertical resonance is 
a crucial factor which reflects a non-local feature of the 
system: its orbital behavior. Indeed and contrary to col- 
lisional fluids, collisionless fluids may develop long range 
correlations which are not captured by a purely local de- 
scription. The distinction between kinematic based and 
spatial mass distribution based i nstabilities in collision - 
less system has been presented bv lPfennigerl (|1996lll998D . 
A fire-hose instability belongs to instabilities depending 
on a strong gradient in the velocity part of the DF, while 
a bar buckling instability belongs to instabilities mainly 
related to the presence of a strong resonance, which is 
determined by the spatial mass distribution. 

In this section, we investigate the role played by the 
shear stress, in particular, the meridional component 
{ T s,rz) which exerts a torque in the vertical direction on 
an imaginary cube of the stellar fluid. 

First, we quantify the buckling amplitude {A\ z ) by 
computing the m = 1 Fourier component in the r — z 
plane of the disk with the major axis of the bar aligned 
to the disk major axis and identify the buckling phase by 
studying the time evolution of A\ z . During the the buck- 
ling phase A\ z reaches a peak value and sometimes goes 
throu gh a second buckling ([Martinez- Valpuesta et al.l 
120061) . In Fig. [71 we present the time evolution of A\ z for 
the three models. A\ z for RCG051A shows a strong peak 
at T ~ 2 Gyr usually considered as the first buckling of 
the bar. It is interesting to note that the bar in this 
model suffers subsequent buckling of smaller amplitudes. 
The onset of the buckling instability can be indicated 
by different physical parameters, e.g., a drop in A2 or 
<r z /a r as mentioned above. For the model RCG051A, 
both Fig. [5] and Fig. [8] indicate that at around 2 Gyr, 
there is a drop in A2 and a z ja r respectively as found in 
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Fig. 8. — Time evolution of the flattening of the velocity ellipsoid 
in galaxy models RCG051A(green solid line), RHG097 (red dashed 
line) and RHG057 (blue dotted line). 
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Fig. 9. — Time evolution of the meridional shear stress nor- 
malized by the initial normal stress along the radial direction in 
RCG051A(green), RHG097 (red) and RHG057 (blue). 




R / R d 

Fig. 10. — Radial variation of the disk midplanc in the bar region 
for the model RHG097. The profiles are drawn at T = (dash- 
dot-dot line), 3.4 (solid black line), 4.0 (red dashed line), 4.4 (green 
dash-dot line) and 5.5 (blue dotted line) Gyrs. 
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Fig. 11. — Edge-on projection of the surface density maps of the 
stellar disk in model RCG051A. Time shown in the panels is in 
Gyr. Buckling instability is evident from the vertical asymmetry 
of the density contours at t = 2 Gyrs. 



previous studies. 

In Fig. [5J we show the time evolution of the merid- 
ional shear stress normalized by the initial normal stress 
defined in Eq. [3] It demonstrates clearly that the merid- 
ional stress develop slowly as the bar evolves and reaches 
its first peak value as the bar enters the buckling phase 
at about 2 Gyr following closely the time evolution of A\ z 
in RCG051A. We establish, here, a new indicator of the 
bar buckling instability that correlates well with other 
indicators mentioned above in a galaxy with a cold stel- 
lar disk undergoing a rapid phase of bar growth (here, 
RCG051A). Let us now examine the other two models 
where the bar growth rate is rather slow in comparison 
to RCG051A. 

In model RHG097, the peak in A\ z (see Fig. [7]) coin- 
cides with that in the meridional stress shown in Fig. |H] 
at around 3.4 Gyr. However, the drop in a z ja r occurs 
noticeably earlier at around 1.5 Gyr, when no buckling 
event was found from the time evolution of A\ z as well as 
from the visual inspection of the surface density maps in 
edge-on projection (see Fig. I12|) . It is evident that a drop 
in o z ja r is not an unambiguous indicator of the buck- 
ling event of strong bars, whereas the meridional stress 
is promising in indicating the onset of buckling instabil- 
ity. 

On the other hand, for a dark matter dominated radi- 
ally hot stellar disk as in model RHG057, the bar grows 
on a much slower rate (see Fig. [5]) and shows no peak 
in A\ z (see Fig. [7]). Also Fig. |8] shows no appreciable 
drop in a z / ay except at around 6 Gyr where a smooth 
decrease in o z ja r is apparent. The meridional stress re- 
mains nearly flat and close to zero for this galaxy model 
which has been evolved for about 12 Gyr during which 
no buckling event was detected. 




Fig. 12.— Same as in Fig. [TT] but for model RHG097. Buckling 
instability occurs at t = 3.4 Gyrs. 




Fig. 13.— Same as in Fig.[TT]but for model RHG057. No vertical 
asymmetry in the density contours detected within t ~ 12 Gyrs. 

As the bar evolves through the buckling phase, the 
disk midplane also responds and exhibits a characteris- 
tic buckling mode. We compute the location of the disk 
midplane using the following formula to follow the buck- 
ling: 



/ zp{r, z)dz 
J p(r, z)dz 



(5) 
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where p(r, z) is the volume density distribution of stars. 
Since the meridional shear stress for the model RHG097 
was comparatively high, the disk midplane was expected 
to show noticeable bending. We discuss here only the 
case of RHG097 and mention briefly the other models. 

Initially, the disk midplane remains flat at z = as 
shown in Fig.EJfor the model RHG097. At t = 3.4 Gyrs, 
the midplane reaches its peak value (z) ~ 130 pc. Sub- 
sequently, the z-amplitude decreases to nearly zero at 
around 5.5 Gyrs restoring the symmetry along the ver- 
tical direction. The buckling modes of the bar in this 
model have characteristic nodes at R ~ Q.6Rd and 
~ 1.8i?d (see Fig. [TU]). Comparing with Fig. [51 we find 
that the location of the second node is close to the coro- 
tation of the bar. The time evolution of the z-amplitudc 
indicates that the buckling instability is a sudden event 
in the galaxy evolution. 

The model RCG051A also showed similar behaviour in 
the ^-amplitude. But the ^-amplitude in model RHG057 
remained close to zero at all times. 

5.1. Meridional tilt of the velocity ellipsoid 

In this section, we discuss the orientation of the stellar 
velocity ellipsoid in the meridional plane (R — z plane) 
of the galactic disk. Since the buckling instability cre- 
ates an asymmetry in the vertical density distribution 
and assuming it preserves reflection symmetry with re- 
spect to the galactic centre, we consider only one half of 
the meridional plane for the computation of the stellar 
velocity ellipsoid. When the bar has formed, we rotate 
it so that its major axis is perpendicular to the line-of- 
sight; in other words, the tilt is calculated in the merid- 
ional plane containing the bar major axis. In edge-on 
projection, the meridional plane would closely resemble 
the surface density maps shown, for example, in Fig. 1121 
In order to understand the spatio-temporal variation of 
the stellar velocity ellipsoid in a model galaxy, we fur- 
ther subdivide the entire meridional plane into several 
small cells each of which contain sufficient number of 
particles for reasonable estimate of the velocity disper- 
sion and the meridional tilt. The cell sizes are fixed at 
AR = QARd in radius and Az = 0.2Rd- The number 
of particles in each of these cells vary over time as they 
are subject to mixing and migration driven by the com- 
bined effie^t_of_aji_evolving bar a nd spiral structures in the 
disk (jMinchev fc Famaevi [2010h . To give an idea of the 
number of particles used, the cell at R = 1.8 — 2ARd 
and z = 0.3 — 0.5Rd contains about 20, 000 particles 
and the innermost cells have about 100000 particles at 
t = 4.8 Gyr for the model RHG097 (for reference see 
Fig. [Ml). 

Fig. [H] depicts the spatio-temporal variation of the 
velocity ellipsoid in the meridional plane of the galaxy 
model RHG097. Initially, the velocity ellipsoid in the 
disk are all aligned with the galactocentric coordinate 
axes and the same holds true for a period of about 3 Gyr 
when the bar has fully developed in the disk, see Fig. [3l 
At 3.4 Gyr, the bar undergoes a sudden buckling instabil- 
ity and the meridional tilt of the velocity ellipsoid reaches 
a peak value as can be seen from the corresponding panel 
in Fig. [TJ] Note that the maximum of the tilt occurs 
at the first node of the buckled bar (see Fig. ITU)) which 
roughly coincides with the edge of the peanut shape in 
this model (see Fig. [12]). In general, higher values of tilt 



angle can be found in the b/p region away from the mi- 
nor axis of the galaxy during the buckling phase. It is 
interesting to notice the spontaneous symmetry breaking 
in the shape distribution of the velocity ellipsoid in the 
meridional plane about the midplane of the galaxy just 
after the peak of the buckling phase. Such asymmetry 
continues to persist for about 1 — 1.5 Gyr since the onset 
of buckling instability, during which the density distribu- 
tion is also asymmetric about the midplane (see, Fig.[T2|). 
After the peak of the buckling phase, the tilt angle of 
the velocity ellipsoid gradually decreases to a low value 
during the subsequent evolution of the galaxy, restoring 
symmetry both in the shape distribution of the velocity 
ellipsoid and the mass density about the midplane. 

Note that, the velocity ellipsoid near the minor axis of 
the galaxy remains nearly aligned with the galactocentric 
coordinate axes before the onset of buckling instability 
and at later times. Although not shown explicitly here, 
the meridional tilt angle of the velocity ellipsoid along the 
minor axis of the galaxy remains zero at all times during 
the galaxy evolution. The meridional tilt of the veloc- 
ity ellipsoid outside the b/p region is nearly unaffected 
by the buckling instability. As shown clearly in Fig. [JJ] 
(see panel at t = 2.4 Gyr), the meridional tilt angle is 
nearly zero for galaxies which host a bar that did not go 
through a buckling instability. From Fig. Q31 it is clear 
that the size of the velocity ellipsoid near the minor axis 
nearly doubles at times when buckling instability is at its 
maximum and their sizes continue to increase. Since the 
semi-major axis of the velocity ellisoid actually measures 
the radial velocity dispersion, it shows clear indication of 
heating in the whole b/p region of the galaxy model. 

Fig. [TBI shows the time evolution of the average merid- 
ional tilt angle (©tilt) of the velocity ellipsoid computed 
in the b/p region for all the models. In both models 
RHG097 and RCG051A, the meridional tilt angle rises 
to a peak value during the buckling instability phase. On 
the other hand, the tilt angle scatters around zero for the 
model RHG057 at all times. From the time evolution, a 
large value of the tilt angle is a characteristic signature of 
the buckling phase that these model galaxies might have 
undergone. In other words, findings of a large value of 
the tilt angle in the b/p region of a galaxy would indicate 
that it might be in the buckling phase or near the vicinity 
of the buckling instability. The subsequent evolution of 
the buckling instability in model RHG097 is particularly 
interesting because of the gradual decrement of the tilt 
angle. It takes about 1 Gyr for the meridional tilt angle 
to fall by half its peak value and can be considered as the 
half-life of the buckling phase (Ttilt,i/2) the galaxy has 
experienced. We find that T t ii t l / 2 ~ 3x T or b, where T or b 
is the orbital time at the disk half-mass radius (for this 
model), which is quite short compared to the galaxy's 
lifetime. This might be the reason for the difficulty in ob- 
serving galaxies in the buckling phase. However, T tilt l / 2 
may depend on various parameters of the galaxy models 
and a thorough search of the parameter space is required 
to find an optimal galaxy model which would show large 
values of tilt angle over long periods of time. The depen- 
dence of Ttiit,i/2 on the dark halo and bulge properties 
will be considered in a future paper. 

5.2. Second moment of DF 
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Fig. 14. — 2D map of the stellar velocity ellipsoid in the meridional (R — z) plane of the disk in model RHG097. Time units are in Gyrs. 
Color bar represents the amplitude of the radial velocity dispersion (oy). The actual value of oy is obtained by multiplying the color indices 
by 893 kms^ 1 . The major and minor axes of the velocity ellipsoid are deter mine d by oy and a z and they are denoted as inscribed crosses. 
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Fig. 15. — Time evolution of the meridional tilt angle of the 
stellar velocity ellipsoid in the boxy bulge region in three models- 
RCG051A (green), RHG097 (red) and RHG057 (blue). 

The calculation of the stellar velocity ellipsoid assumes 
that the moment integrals of the DF exist and returns 
a finite value. However, the validity of such an assump- 
tion is questionable, especially when the stellar system 
is undergoing an unstable phase e.g., buckling instabil- 
ity in the present case. The resonant parts of the phase 
space during such an instability can develop bi-modal 
and/or particle distribution with long tail for which the 
very notion of first or second moment of the DF is math- 
ematically no longer meaningful. Bimodal velocity dis- 
tribution has been observed for late-type stars in the so- 
lar neighborhood by Hipparcos and numerical models of 
disk response to a bar is shown to have reproduc ed many 
such features in the local velocity dis tribution (jDehnenl 
120001 : |Fu3 I200T1 : IMinchev et a!ll20Tnl) . In Fig. QH we 
show the normalized histograms for radial (v r ) and verti- 
cal (v z ) velocities in the model RHG097 during and after 
the buckling instability nearly disappeared. We picked 
up three different regions in the meridional plane and 
histograms, in the two regions (R = 0.6, z = 0. and 
R = 0.6, z = 0.2) where the meridional tilt was max- 
imum, are fairly well represented by a single Gaussian 
DF with different variances. In the region close to the 
minor axis of the galaxy i.e., R = 0.2, z = 0.2, the radial 
velocity histograms needed two Gaussian DFs: one with 
cold component with dispersion ~ 30 kms -1 and one 
with a hot component with a dispersion ~ 102 kms . 
A close inspection of Fig. [14] indicates that the stars are 
heated strongly in that region as it is clear from the size 
of the ellipsoids. At the time of buckling, the size of the 
ellipsoid nearly doubles indicating an increase in the ve- 
locity dispersion by a factor of two. In any case, in all 
the regions examined, we have a unimodal DF to repre- 
sent the stars in the meridional plane and they show well 
behaved first and second moments. 

6. DISCUSSION AND CONCLUSIONS 

The buckling instability is one of the routes through 
which an initially axisymmetric stellar disk would form a 
boxy/peanut bulge such as present in our own Galaxy. In 
order to understand the formation of such a boxy /peanut 
morphology, it is important to have further insight on the 
buckling instability. How and when would a bar go buck- 
ling unstable? How many buckling events has a present 
day galaxy experienced? How does it depend on the 
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Fig. 16. — Velocity histograms of stars at 3 different regions (as 
indicated in the panels) in the meridional plane of the galaxy model 
RHG097. Blue lines indicate vertical velocity and red and broader 
ones indicate radial. Only the radial velocity histograms arc fitted 
with Gaussians, just for illustration. 

dark matter fraction in galaxies? There are several is- 
sues needed to be addressed in order to grasp this phe- 
nomenon. The current paper addresses one such issue on 
the onset of a buckling event in a disk galaxy. It is shown 
that there is a connection between the onset of the buck- 
ling instability and the shear stress of stars. We see that 
the shear stress reaches its peak value during the buckling 
phase and then decreases gradually. The development of 
a shear stress in the stars is a result of collective pro- 
cess in the disk. If these stars are also trapped in the 
vertical IL Rs, this would eventually lead to the buc kling 
instability (jPfenniger k FriedlilH99ll lQuillenll2002D . 

From our study, it is clear that a bar that grows very 
slowly, on a several Gyr time scale, does neither develop 
any appreciable shear stress nor go through any buck- 
ling instability. On the other hand, bars that grow very 
rapidly such as in RCG051A, develop both, shear stress 
and buckling instability. In no cases that we have stud- 
ied does a shear stress develop in the bar and not go 



10 



Saha et al. 



through the buckling instability. One emerging scenario 
is that the development of shear stress is related to the 
rate at which a bar grows i.e., the rate at which a bar 
strength grows through transport of angular momentum 
outward. In galaxy models with higher values of Toomrc 
Q, the growth of bar strength is rather slow leading to 
weak bar and insignificant amount of shear stress. 

Several important kinematical changes occur in the 
galaxy during and after the episode of the buckling in- 
stability. From the Fig. HU it is clear that the stars are 
heated in the b/p region, especially close to the minor 
axis of the galaxy, by a factor of ~ 2 durin g the buck- 
ling in stability. This was shown previously bv lSaha et al.l 
•2010: for model RHG097 and others. Another inter- 
esting aspect is the kinematical changes in the vertical 
structure of the galaxy. Note, the stellar disk is isother- 
mal initially (see the first panel at t = 1.1 Gyr in Fig. [T4")) . 
In the after-episode of the buckling phase, there is a 
clear distinction in the velocity dispersion above and be- 
low the disk midplane indicating spontaneous breaking 
of isothermal structure in the b/p region. Such non- 
isothermal vertical structure in the b/p region is persis- 
tent long after the buckling phase. We will address this 
issue in more detail in a future paper. 

Our main conclusions from this work are as follows: 
1. We show that the meridional tilt of the stellar veloc- 
ity ellipsoid is a better indicator compared to a drop in 
the bar amplitude or cr z /a r for the onset of the buck- 
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